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The scaling equation of state of the three-dimensional O(N) universality 
class: N > 4. 
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We determine the critical equation of state of the three-dimensional O(N) universality class, for N = 4, 5, 
6, 32, 64. The N — 4 is relevant for the chiral phase transition in QCD with two flavors, the N = 5 model 
is relevant for the SO(5) theory of high-T c superconductivity, while the N = 6 model is relevant for the chiral 
phase transition in two-color QCD with two flavors. We first consider the small-field expansion of the effective 
potential (Helmholtz free energy). Then, we apply a systematic approximation scheme based on polynomial 
parametric representations that are valid in the whole critical regime, satisfy the correct analytic properties 
(Griffiths' analyticity) , take into account the Goldstone singularities at the coexistence curve, and match the 
small-field expansion of the effective potential. From the approximate representations of the equation of state, we 
obtain estimates of universal amplitude ratios. We also compare our approximate solutions with those obtained 
in the large- TV expansion, up to order 1/iV, finding good agreement for N > 32. 



1. Introduction 

In the theory of critical phenomena, continuous 
phase transitions can be classified into universal- 
ity classes determined only by a few properties 
characterizing the system, such as the space di- 
mensionality, the range of interaction, the num- 
ber of components of the order parameter, the 
symmetry and the symmetry-breaking pattern 
PP. Renormalization-group theory predicts that 
critical exponents, universal amplitude ratios and 
scaling functions are the same for all systems 
belonging to a given universality class. The 
O(N) universality classes are among the most 
important ones. They are characterized by an 
TV-component order parameter and a symmetry 
group 0(N) which is spontaneously broken to a 
subgroup 0(N — 1) in the low-temperature phase. 
For a recent review on this subject, see ref. [2]. 

Here we study the O(N) three-dimensional 
universality class, with N > 4. The three- 



dimensional 0(4) model is relevant for the finite- 
temperature behavior of QCD with two light- 
quark flavors [3]. The 3-D 0(5) model is rel- 
evant for the so-called SO (5) theory of high-T c 
superconductivity jjjj. According to universality 
arguments, the N — 6 model should describe the 
chiral phase transition in two-color QCD with two 
flavors 

The equation of state is a relation between 
the magnetization M, the reduced temperature 
t = (T — T c )/T c and the external magnetic field 
H. Near the critical point it has the scaling 
form [U 



H = {B c )- 5 MM°- L f(x) 



6-1 . 



(1) 



where M = \M\, x = B l /l 3 tM^ 1 / is a scaling 
variable, f(x) is a universal scaling function, fixed 
by the normalizations /(0) = 1, /(— 1) = 0, B c 
and B are the non-universal magnetization am- 
plitudes at the critical isotherm and at the coex- 
istence curve: 
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M = B C HH 1/S ~ 



t = 0, 



(2) 
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M = B(-t) 







t<0, H^O. 



(3) 



We introduce two approximation schemes: 



The equation of state can also be written in the 
form: 



M 

H = ab—t 0S F(z), zsbMt~ p , 
M 



(4) 



where the non-universal constants a and b fix the 
normalization on F(z) such that for z — > 



r 2n _ z 2n-l 



;i>3 



(2n-l)! 



(5) 



Inverting the scaling equation of state can be 
written as: 



M = B c HH x ' 5 - x E{y), 



(6) 



where y = [B / Be) 1 * '^tH~ x /^ s > is a scaling vari- 
able. All the functions f(x), F(z) and E(y) are 
universal. 

2. Approximate equation of state 

The parametric representation 
M = m R /3 m(9) 



t = R(l-8 2 ) 
H = h {) R f)S h{9) 



(7) 



(mo and ho are normalization constants) imple- 
ments the known analytical and scaling proper- 
ties of the equation of state. R is a nonnega- 
tive variable which measures the distance from 
the critical point. The functions h(0) and m(9) 
are odd and are conventionally normalized so 
that h{6) = 9 + 0(9 3 ) and m(6) = 6 + 0((9 3 ), for 
9^0. As can be seen from JJJ, the line 9 — 
corresponds to the high-temperature phase, while 
on the critical isotherm 9 = 1, The coexistence 
curve is given by 9 = 9q, where 8o is the first 
positive zero of h(9). Near x = — 1 (coexistence 
curve), f(x) ~ Cf(x + l) 2 . This can be satisfied 
if h{9) - (6»-6» ) 2 , for 0-> 6» . 

We start from the high-temperature, small 
magnetization, expansion z — > of eq. JSJ. 

Then, using the representation J7J), we perform 
an analytic continuation to the low-temperature 
phase. 



A : < 



B : { 



m(9) = i 

h{9) = i 
( m{8) = 

h{8) = 



1 + E 



»=i 

fl2\ 2 



(8) 



For n = the two schemes are the same. In 
both cases, the coefficients Cj and # are deter- 
mined by imposing that the equation of state, for 
z — > or, equivalently, 9 — > 0, reproduces the 
expansion JSJ. We have considered both schemes 
for n = 0, 1. We use field-theoretical estimates 
of r6 and rg, obtained by analyzing the pertur- 
bative series |6I7) . For the critical exponents we 
use the Monte Carlo estimates of [S] for N = 4 
and field-theoretical estimates [7] for the other 
models. Consistence of the whole computation 
requires the coefficients Cj to be small. Moreover, 
the Jacobian of Q must not vanish in the interval 

[0M 

3. Results 

We show in figure 2] taken from the scal- 
ing function f[x) for the 0(4) model, as obtained 
with the n — 0, n — 1 A and n = 1 B schemes. 
We also show a comparison with the Monte Carlo 
result of ref. ^01- For N > 4, scheme B did not 
work, failing to satisfy consistency conditions (see 
end of section In figure [3 taken from [7], we 
show the scaling function E(y) for the 0(5) and 
0(6) models; n — 1 refers to scheme A. For the 
0(6) model, we also show a comparison with the 
Monte Carlo result of ref. JT] . 

From the approximate scaling equation of 
state, it is possible to obtain several amplitude 
ratios. We report in tabled two of the most im- 
portant ones: the ratio of the amplitudes of the 
specific heat U and R x = C + B 6 ' 1 / (B c ) 5 , where 
C + is the amplitude of the susceptibility in the 
high-temperature phase We also show the 
same quantities for the N = 32, 64 models, along 
with a comparison with the large-TV expansion. 
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Table 1 

Universal amplitude ratios. 



N U 



4 


1.91(10) 


1.12(11) 


■5 


2.2(2) 


1.2(1) 


6 


2.5(2) 


1.15(9) 


32 


1.5(5) 


0.94(1) 




1.47* 


0.940* 


64 


1.5(4) 


0.968(4) 




1.47* 


0.9701* 



* Large- TV 



Figure 1. The scaling function f(x) for the 0(4) 
model. A comparison with the Monte Carlo result 
of ref. ^Oj is shown. From [3J. 




Figure 2. The scaling function E(y), for models 
N = 5 (up) and N = 6 (down). For the 0(6) 
model, we also show a comparison with the Monte 
Carlo result of ref. ^J. From [7J. 



In the N — > oo limit, the n — scheme becomes 
exact ^21; f° r the 0(32) and 0(64) models we 
find a good agreement with the large- N results. 

In conclusion, we note that the results obtained 
in the n = 1 schemes are close to those in the 
n = scheme, thus supporting the effectiveness 
of the approximation scheme. For a more de- 
tailed discussion, as well as for other results, see 
refs. [9TT| . 
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